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In this talk I'm going to summarize work done with Larry Yaffe jl] on applying e 
expansion techniques to the electroweak phase transition. As I shall review, analysis 
based on the perturbative expansion of the effective potential, which has already been 
discussed in this conference, is only valid (in the minimal standard model) when the 
zero-temperature Higgs mass is small compared to the zero-temperature W mass. The 
goal of my talk will be to discuss what can be done in the experimentally more realistic 
case when this condition fails. I shall discuss how one might compute various parameters 
of the transitions such as the correlation length at the critical temperature, the latent 
heat of the transition, the bubble nucleation rate, and the baryon violation rate after 
the transition. Note that these are all quasi-equilibrium quantities. In particular, I 
shall modestly not attempt to tackle the much harder problems directly associated 
with the generation of baryon asymmetry on the expanding bubble walls. 

Throughout this talk I will work in a simple toy model: the minimal standard 
model with a single Higgs doublet. I do this because the model has relatively few 
unknown parameters, and so it is relatively easy to analyze and results are not obscured 
by complicated parameter dependence. I call it a toy model because I don't care to 
step into the debate of whether the minimal model has enough CP violation |2], [|. For 
simplicity, I will also be ignoring the Weinberg angle and focusing on a pure SU(2) 
Higgs theory. 

Why life is simple when m(Higgs) T=0 <C M(W) T=0 

* Talk presented at the NATO Advanced Workshop on Electroweak Physics and the Early Universe: 
Sintra, Portugal, 1994. 



Consider the classical Mexican-hat potential for the Higgs: 

V ~ -/i 2 2 + A0 4 . (1) 

And now consider the free energy of the system at finite temperature in the background 
of a Higgs field 0. One contribution will be from the classical energy Vq above for that 
0. But there will also be contributions from the free energy of the real particles present 
in the plasma. Focus on the contribution from W bosons. This contribution (at one- 
loop) is nothing more than a formula you can look up in your old graduate statistical 
mechanics book; it is the free energy of an ideal Bose gas: 



AF ~ T / d 3 k In (l - e~ pEk ) , (2) 

where the energy of a relativistic particle is 



E k = VP + M 2 ,- yjk* + g*<t>*. (3) 

Note that the W gets its mass from the background Higgs field, so that M\y ~ 9<P- 
Thus AF above depends on the background value of 0. At high temperatures, one 
can expand AF in powers of My?/T, and one gets a high-temperature expansion of the 
form 

AF ~ #T 4 + #M^T 2 - #M^T + • ■ • , (4) 

where # indicates numerical constants that I'm not going to bother writing down. 
Putting the classical potential V together with the W contribution to AF givesjj] 

F ~ V + AF ~ const + (-/j 2 + g 2 T 2 )(f) 2 - g 3 (f) 3 T + A0 4 + ■ • ■ (5) 

(where I'm no longer going to bother even putting in the # signs). The g 2 T 2 <p 2 term 
comes from the M 2 T 2 term of the high-temperature expansion of AF. This term is 
responsible for the phase transition by turning the curvature of the potential at the 
origin from negative at low temperature to positive at high temperature. The g 3 <p 3 T 
term, which comes from the M 3 T 2 term, is responsible for making the phase transition 
first-order as shown in fig. [l|. Specifically, there exists a temperature T c at which 
there are two degenerate minima of the potential, and the expectation of changes 
discontinuously as one cools through the transition. I shall ignore the generally small 
effects of subsequent terms in the high-temperature expansion (Q). 

Now consider the T = T c curve in fig. |l|. The shape comes from the interplay of 
the 2 , 3 and 4 terms in the free energy (|5|): the 2 term turns it up at the origin, the 
(f) 3 term turns it down again, and the 4 term turns it up at large <fi. In particular, in 
the region of <j) near the symmetry-breaking minimum (where (f) ^ 0) or near the hump 
in the potential, these three terms must be the same order of magnitude, so that 

(-/. 2 + ^TV~^ 3 T~A0 4 . (6) 

From the last relation in particular, one finds that this region of is characterized by 

q 3 
0~yT. (7) 

Now consider the order of the loop expansion parameter if we perturbatively an- 
alyze the physics of the phase transition. Each loop costs a factor of g 2 . However, 



t Of course, there are also scalar and other contributions to AF. I'm just focusing on the W 
contributions here for the sake of brevity and pedagogy. 




Figure 1: The form of the free energy, as a function of 0, for different temperatures. 



at high temperatures, each loop is also associated with a factor of the temperature 
T, so that the cost of a loop is actually g 2 T. (One way to see this factor of T is to 
remember that, at finite temperature, Euclidean time becomes periodic with period 
(3 = 1/T. As a result, zero-temperature loop momentum integrals J d A p are replaced 
by finite Fourier sums T J2 P0 1 d 3 p over the frequencies, and a factor of T then appears 
explicitly.) Finally, by dimensional analysis, g 2 T must be divided by the mass scale 
Mw of the problem to determine the final, dimensionless cost of adding a loop. The 
loop expansion parameter (after resummation) is then 
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where I have used Mw ~ g<t> and (0). A and g 2 , however, are related to the zero- 
temperature Higgs and W masses, so that the loop expansion parameter can be rewrit- 
ten as 

X_ m 2 (Higgs) T =o , v 

g 2 ~ M 2 (W) T=0 ' l ' 

So we see that the perturbative loop expansion is only valid when the zero-temperature 
Higgs mass is small compared to the zero-temperature W mass. All of this power 
counting is reviewed in ref. [01. 



Why life is not simple 



One of the constraints of electroweak baryogenesis is that there must not be any 
significant amount of baryon number violation after the completion of the phase tran- 
sition. Otherwise, any baryon asymmetry generated during the transition would be 
washed away to zero. The rate of baryon number violation is exponentially sensitive 
to the sphaleron mass: 



rate ~ e -^ s P haleron ) ~ e 



#M w /g 2 T ^ e -#g 2 /X 



(10) 



The second equality above comes because the sphaleron mass is of order Mw/«„, and 
the final equality because the exponent My^j g 2 T is (not coincidentally) nothing other 
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Figure 2: The Boltzmann exponent for baryon number violation vs. the zero- 
temperature Higgs mass. 



than the inverse of the loop expansion parameter (|8|), which we previously determined 
for the broken-symmetry phase. You can now see that the rate will be small after the 
transition only if the zero-temperature Higgs mass (and hence A) is small. The inverse 
relationship between the rate exponent and the zero-temperature Higgs mass is shown 
schematically in fig. ^. Requiring that the rate be small compared to the expansion 
rate of the universe puts a lower bound on /^(sphaleron) and hence an upper bound 
on the Higgs mass. A careful analysis, based on the one-loop effective potential, yields 
an upper bound on the Higgs mass of 35-40 GeV due to this constraint . This seems 
to rule out electroweak baryogenesis in the minimal standard model. 

But now it's time to wonder whether a one-loop analysis can be trusted. In par- 
ticular, is a Higgs mass of 35 GeV in fact small compared to the W mass of 80 GeV? 
That depends on all the factors of 2 left out of the order-of-magnitude analysis of loop 
corrections that I presented above. To settle this issue, the two-loop potential was com- 
puted in ref. 0] (and the most dominant pieces independently in ref. |J). Fig. |3| shows 
the difference between the one-loop and two-loop potentials at the critical temperature. 
Whereas the VEV changed only by about 20%J^] the height of the hump changed by 
a factor of three! The moral is that, for some quantities, (35 GeV)/ (80 GeV) is not a 
small number. And so it is not clear, in particular, whether higher-loop corrections to 
the baryon number violation rate will be small. 

But what can we do when mean field theory (i.e. perturbation theory) fails — 
that is, when fluctuations around the mean field <fi are large and cannot be treated 
perturbatively. To phrase it another way, what do we do when the correlation length 
1/Mw becomes very large at T c so that the loop expansion parameter (||) is large? 
Condensed matter physicists have long dealt with this problem in the context of second- 
order phase transitions, where the correlation length is infinite at T c . Our strategy is 
to borrow one of their techniques — the e expansion — which has been very successful for 
some systems.^] 

* Beyond leading order the VEV is an unphysical and gauge-dependent quantity. The result shown 
here is for Landau gauge. 

§ The idea of applying the e expansion to the electroweak phase transition has been suggested 
previously by Gleisser and Kolb |m], March-Russel ||, and probably by others that I'm not familiar 

with. 
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Figure 3: The (resummed) two-loop vs. one-loop free energy, as functions of 0, for 
m(Higgs)x , =o — 35 GeV at the critical temperature (as determined individually for each 
potential). 



The e expansion 

The e expansion consists of the following steps. 

Step 1. Recall that equilibrium quantities can generally be expressed in terms of 
the partition function Z = tre _/3H . The exponential e~^ H looks just like an evolution 
operator e tHt when the time t — if3 is imaginary. Just like the evolution operator, Z has 
a path-integral representation, and the only difference with standard zero-temperature 
path-integrals is that time only extends for Euclidean time f3: 



[D(f>] exp 



drLi 



IT 



Also, the trace is implemented by integrating over all configurations where the initial 
state is the same as the final state, which translates into the periodic boundary con- 
dition 0(0, x) = (f>(/3, x) on the path integral. (Fermions have anti-periodic boundary 
conditions.) 

Step 2. Since the Euclidean time dimension only has extent j3, it disappears in 
the high temperature limit (3 — > 0. That is, if one is studying (equilibrium) physics at 
distance scales I large compared to 1/T, then the Euclidean time dimension decouples 
and the four-dimensional theory is reduced to an effective three-dimensional theory of 
the static modes (po — 0). The effects of fluctuations (pq ^ 0) in the Euclidean time 
direction decouple, as powers of (3/1, except for renormalizations of relevant operators. 
For example, the three dimensional mass and couplings are related to the original four- 
dimensional ones by 



m. 



fh 



m 4 + g T + higher-order . 
g1(T) + higher-order . 



(12) 
(13) 



The g 2 T 2 term in the mass is just the term I discussed earlier that turns the curvature 
of the potential at = from concave down at low temperature to concave up at 
high temperature, and so drives the phase transition. In four-dimensional language, 
the phase transition is achieved by varying T. In three-dimensional language, that 
corresponds to varying the mass m\ from negative values up through positive values. 




Figure 4: The classical potential for scalar theory, as a function of the mass of the 
effective 3 dimensional theory. 



Step 3. As we shall see, the three-dimensional theory is difficult to directly solve 
in a well-defined systematic fashion. The trick of the e expansion is to solve instead 
a related class of theories by generalizing the 3 spatial dimensions to 4 — e spatial 
dimensions, replacing the action in the usual way by 



S(d=3) -> / d 



{D<f)) 2 + m 2 <f) 2 + //A0 4 + 



(14) 



(Note that I am defining my couplings, such as A, to be dimensionless.) 

Step 4- For e<l, the theory turns out to be solvable perturbatively in e by using 

renormalizat ion-group (RG) improved perturbation theory. I'll discuss this is more 

detail in a moment. 

Step 5. After computing to a given order in e, return to the world of three spatial 

dimensions by taking e — > 1. 

Review of a success story: pure scalar theory 

Forget about gauge theories for a moment and consider a theory of a single, real 
scalar field: 



U 



(d(j>) 2 + m 2 2 + //A</> 4 



(15) 



The phase transition of the tree-level potential is depicted in Fig. |], which shows a 2nd- 
order transition. The 1-loop potential (appropriately resummed) for this model turns 
out to predict a lst-order phase transition, but the loop expansion parameter turns 
out to be order 1 and so mean field theory results cannot be trusted.^] In fact, this 
scalar theory goes by a familiar name near its critical point: it is the Ising model, and 
the Ising model is known to have a 2nd-order phase transition. (More accurately, the 
scalar theory is in the same universality class as the Ising model. The field <p roughly 
corresponds to the Ising spin averaged over a volume of size l/// 3 .) 

In a 2nd-order phase transition, the correlation length becomes infinite. To study 
such a transition, one should therefore examine the system on larger and larger distance 



" For a review, see section II. A of ref. 
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Figure 5: The renormalization group flow, into the infrared, of the pure scalar theory. 



scales and try to determine its scaling behavior. But, as particle physicists, we're all 
quite familiar with how to change scale in 4 — e dimensions. We use the renormalization 
group: 

fjtdpX = -eA + cX 2 + 0{\ 2 ) . (16) 

If one set e to zero, this would just be the usual four- dimensional renormalization group 
for a pure scalar theory, with the one- loop /^-functions depicted explicitly as cX 2 . c is 
a numerical constant that is easily computed. Slightly away from four dimensions, the 
additional term — eA on the right-hand side simply represents the classical scaling of 
the interaction A0 4 due to its dimension. For example, at tree level, a change in /i by a 
factor of 2 can be compensated in //A0 by changing A by a factor of 2~ e . This trivial 
scaling is the sole content of the — eA term in fll6"|). It is the other terms which contain 
the physics of integrating out degrees of freedom as the scale \x is changed. 

The flow of the RG equation (|HJ) is shown in Fig. |5|. The arrows depict flow into 
the infrared (/i — > 0). At small A, the — eA term dominates and drives A bigger; at 
large A, the +cA 2 term dominates and drives it smaller. There is an infrared-stable 
fixed point when the two terms balance at A* = e/c. The long-distance behavior of the 
system will be described by this fixed point. 

When e is small, the fixed point A* = e/c will be at small coupling, and therefore 
our neglect of two-loop and higher-order contributions to the renormalization group 
in ( |T5] ) will be justified. More generally, the loop expansion in A at the fixed point 
is equivalent to an expansion in e. We know, however, that the loop expansion does 
not converge but is asymptotic, blowing up at order A™ where n ~ 1/A*. That means 
that the expansion in e will also be asymptotic, blowing up at order e n where n ~ 1/e. 
If we're lucky, n ~ 1/e will mean 3 or 4 when e — > 1 and the first few terms of the 
expansion will be useful in three dimensions, giving a reasonable approximation for 
quantities of interest. If we're unlucky, n ~ 1/e will mean zero when e — *■ 1, and the 
asymptotic series will start blowing up immediately after the first term and be useless 
in three dimensions. 

Are we lucky or unlucky? Let's look at some examples. The first is the suscepti- 
bility exponent 7 defined by x ~ \T — T c | 7 where x measures the response of (4>) to an 
external source J<p. The result is 

7 = 1 + 0.167e + 0.077e 2 - 0.049e 3 + • ■ • , (17) 

which, for e = 1, gives 1.195 + • ■ -. This compares favorably to the result 1.2405 + 0.0015 
which comes from more sophisticated techniques (the e expansion combined with Pade 
resummation, or renormalization group studies directly in three dimensions). It is also 
consistent with lattice results, which have larger error bars. It is also consistent with 
experiments, which have yet larger error bars. I should also give a much more marginal 
example, which is the correlation exponent r] defined by (0(r)0(O)) ~ r 2 ~ d -v as r _► q^ 
r) is simply twice the anomalous dimension of <j) and is given by 

77 = 0.0185e 2 + 0.0187e 3 - 0.0083e 4 + 0.0359e 5 + ■ • • . (18) 




Figure 6: The renormalization group flow, into the infrared, of electroweak theory. 



If one adds the first three terms for e = 1, and ignores the fourth where the series is 
clearly blowing up, one gets 0.0289. The actual result is believed to be about 0.035. So 
the naive application of the e expansion gives an approximation to within about 20%. 

Electroweak theory 

The one-loop renormalization group equations for the couplings of electroweak 
theory are similar to the case of scalar theory: 






-eg 2 + p g\ 

-eg 2 + (ag A + bg 2 X + c\ 2 ) 



(19) 
(20) 



The scalar /3-function is a little more complicated now because there are terms involving 
g 2 as well as A. 

When one computes the numerical coefficients in these equations and plots the 
renormalization group flow, one gets Fig. |(| The g 2 = line is similar to our previous 
result for pure scalar theory. Howvere, because of the flow in g 2 , the fixed point is 
no longer infrared stable. As observed by Ginsparg ||, there are in fact no infrared 
stable fixed points seen in the e expansion. This suggests that the phase transition 
might always be lst-order. Another, extremely hand-waving, suggestion of this comes 
from the observation that all theories flow into the region A < and, at tree level, 
A < would be an unstable theory. This suggests that there is something suspect 
about flowing to arbitrarily large distances, as one would desire if the transition were 
2nd-order. 

The fact that the transition is lst-order can be made rigorous, for small e, by 
remembering that we know how to solve the theory when A <C g 2 . When X ^ g 2 , 
straightforward perturbation theory is under control and predicts a lst-order transition. 
So now consider a theory initially defined by the couplings g\ and Ai as depicted in 
Fig. |6]. The renormalization group allows us to flow to larger distance scales and fins a 
completely equivalent theory where A <C g 2 , also depicted in the figure. This equivalent 
theory can then be solved with straightforward perturbation theory and predicts a lst- 
order transition. In summary, our method is as follows. 

Step 1. Start with initial couplings gf, Ai at scale fii ~ T (which is the scale where 
the Euclidean time fluctuations first decouple). 

Step 2. Use the renormalization group to flow to an equivalent theory with A -C g 2 . 
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Step 3. We can then perturbatively compute the effective potential or whatever 
other quantities are of interest. 

Step 4- Expand all result in e. 

This procedure is closely related to the work of Rudnick [K| and Chen, Lubensky 
and Nelson Q| 

In the first step, we started with initial couplings g\ and Ai. We know what the 
initial couplings should be in the three dimensional theory; they are perturbatively 
related to the couplings in the original 3+1 dimensional theory as in (JIB]). But in the e 
expansion we have generalized the 3-dimensional theory to a 4— e dimensional theory. 
What couplings g\ and Ai in the 4— e dimensional theory most naturally extrapolate 
to those of the 3-dimensional theory as e — > 1? 

This question is easily answered if one observes that the one-loop renormalization 
group equations are independent of e if appropriately rescaled. Specifically, take 

g 2 ^eg 2 , A - eA , fi - JT 1 ^ . (21) 

The one-loop RG equations (|20|) then become 



fxd.f = -g 2 + [3 g\ (22) 

fid^X = -g 2 + (af + bg 2 \ + c\ 2 ) (23) 

and are independent of e. Thus, if I had simply plotted g 2 /e vs. A/e in Fig. |6| rather than 
g 2 vs. A, the trajectories would have been independent of e. The natural generalization 
of any particular theory from 3 to 4— e dimensions is then clear: just keep the same 
place on the same trajectory. So the relationship between the couplings is 

(^ ) A 4 _ e ) = (e( ? 3 2 ,eA 3 ). (24) 

Note that, even though the trajectories are independent of e once the couplings are 
scaled by e, the rate at which those trajectories are traversed is exponentially sensitive 



to e as a result of the scaling of /x in fl2"T|). In particular, the trajectories are traversed 
exponentially slowly as e — > 0. 

Example 

I'll now very roughly and schematically outline a sample calculation. Consider the 
classical Higgs potential 

V d ~ m 2 (j) 2 + //A0 4 . (25) 

Suppose I've used the renormalization group to run to A <C g 2 . In particular, it's 
computationally convenient to run to A = 0. I can now compute the 1-loop potential, 
and the result for e — > is the familiar Coleman- Weinberg potential: 



01-loop) _ m 2^2 + ^4^4 ln 191 (26) 




The dependence of this potential on the effective mass m gives a lst-order phase tran- 
sition, as shown in Fig. 0. To study the system at the critical temperature, I should 
find the critical value of m where the two minima are degenerate. I can then compute, 

" There is a long story about why the e-expansion prediction that the phase transition is always 
lst-order has, historically, been believed to be incorrect in the U(l) case, where the theory models the 
phase transition of superconductors. For a discussion of problems with and loopholes to this belief, 
see section I.E of ref. H . 




Figure 7: The Coleman- Weinberg potential as a function of m. 



for example, the scalar correlation length £ in the asymmetric phase at the critical 
temperature, which at this order is related to the curvature of the potential. One finds 



£ ~ 



1 



1 



V"( 



fi'V 



(27) 



[The dependence l/g 2 fi 2 is easy to obtain from the form of (p6|) by simple scaling 
arguments.] Now recall from the previous section that g 2 is 0(e) and that the amount 
of running required to run to A = is exponentially long as e — > 0, so that 

(28) 
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(29) 



where # is computed from the renormalization group equations. If one works to higher 
orders in e, one finds that the prefactor of the exponential is a series in e: 



e 2 T 2^I e #A (1 + #e + #e 2 + ...). 



(30) 



When I later refer to orders of the e expansion, I shall be referring to the expansion of 
the prefactor. So "leading-order" will refer to a calculation that includes the first term 
of the prefactor, "next-to-leading order" includes the second term, and so forth. 

Tests 

Now that I've described how to go about computing things with the e expansion, it 
would be nice to find some cases where we can check to see if the e expansion is working. 
One case where we already know how to compute results is in the perturbative regime 
where the initial Xi/gf is C 1. As a first check, we can compute the results of the 
e expansion in the same limit. Table |l] shows the ratio of e expansion results to the 
correct answer in the limit \\/ g\ — > 0. Note that leading-order in e is generally within 
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observable ratio 


LO 


NLO 


asymmetric phase correlation length 


1.14 


0.94 


symmetric phase correlation length 


1.62 


0.92 


latent heat 


0.77 


1.04 


surface tension of a domain wall 


0.60 


0.98 


AF(T ) 


0.24 


0.56 



Table 1: The ratio of the e-expansion results, computing prefactors through leading or- 
der (LO) and next-to-leading order (NLO) in e, to the corresponding three-dimensional 
result when Ai -C q\. 

roughly 50% and next-to-leading order within roughly 5-10%. (The exception is the 
last quantity, which is the difference of free energies of the two ground states at the 
temperature where the = ground state is just loosing its meta-stability.) 

We have also tested the e expansion against large N results where N is the number 
of scalar fields in the theory. I didn't have time to discuss this in Sintra. The result is 
that the e expansion works qualitatively but fails quantitatively at leading-order in e. 
However, the e-expansion self-diagnoses its own quantitative failure because a compu- 
tation of next-to-leading order corrections in e shows that they are large compared to 
the leading-order results. 

The perturbative test discussed above is only of limited interest because, after all, 
we already knew how to solve the theory in the perturbative regime. The large N 
tests are of limited interest because we want to know what happens when N is small. 
(Indeed, the large N theory is very different. For example, it is not asymptotically 
free.) What we really want is a test of the e expansion when Ai > gf. A natural thing 
to do is to check whether the e expansion diagnoses its own failure — that is, whether 
next-to-leading order corrections are small or large compared to leading-order results. 
Computing results at next-to-leading order in this regime is involved, and we only had 
the stamina to make one test. Our test quantity is related to the specific heat. More 
specifically, for technical reasons, it is 



(asymmetric phase correlation length) 2 x (latent heat) 



(31) 



The calculation to next-to-leading order requires the 2-loop renormalization group and 
a 2-loop computation of the potential.^] Fig. |8| shows the results. 

From this figure, we conclude that the e expansion may be quantitatively useful for 
mh below 150 GeV or so, where the correction is within ±30%. The expansion may be 
useful qualitatively, but probably not quantitatively, for larger Higgs mass. You should 
be cautioned, however, that we have only tested one quantity in this way, and perhaps 
it's better behaved than most. 

Qualitative results for B violation 

The original motivation for this study was to examine the bounds on electroweak 
baryogenesis, and so I should say a few words about the implications of the e expansion 
for the rate of B violation at the completion of the phase transition. In particular, one 
can compare the results of renormalization-group improved perturbation theory in 4— e 



** If we had just computed the latent heat by itself, we would have needed the three-loop renormal- 
ization group, which is not currently known. See ref. (T| for details. 
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Figure 8: The relative size of the next-to- leading order correction compared to the 
leading-order result, as a function of zero-temperature Higgs mass. 



dimensions to the results one would get with naive, unimproved perturbation theory 
in 4— e dimensions. At one loop, we find that RG improvement predicts a stronger 
phase transition than the unimproved one- loop result, consistent with the nature of 
the 2-loop correction of Fig. [3] discussed earlier. One might assume that a stronger 
phase transition implies a smaller rate of B violation after the transition. Instead, we 
find a larger rate of B violation in 4— e dimensions. This conclusion is, as always, only 
rigorous for e < 1 and may or may not hold qualitatively when e = 1. (Ideally, one 
should compute next-to-leading order corrections and see if they are under control.) 
However, our result offers an important warning. If one studies the transition on a 
lattice, and finds a transition much stronger than expected by perturbation theory, 
that does not necessarily mean that baryogenesis is more viable in the model being 
studied. 

I should conclude with a final caveat, however. The calculation of the B violation 
rate in the e expansion is beset by additional technical perils arising from the intrinsi- 
cally 3-dimensional nature of the sphaleron. And so the extrapolation of e — > 1 may be 
even more suspect than usual. See ref. QTJ] for details. 
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